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Abstract 

The influence of dark matter particle decay on the baryon-to-photon ratio has been studied for different cosmo¬ 
logical epochs. We consider different parameter values of dark matter particles such as mass, lifetime, the relative 
fraction of dark matter particles. It is shown that the modern value of the dark matter density Ocdm = 0.26 is 
enough to lead to variation of the baryon-to-photon ratio up to A 77/77 ~ 0.01 4-1 for decays of the particles with 
masses 10 GeV 4 -1 TeV. However, such processes can also be accompanied by emergence of an excessive gamma 
ray flux. The observational data on the diffuse gamma ray background are used to making constraints on the 
dark matter decay models and on the maximum possible variation of the baryon-to-photon ratio A 77/77 < 1 CF 5 . 
Detection of such variation of the baryon density in future cosmological experiments can serve as a powerful 
means of studying properties of dark matter particles. 
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1. INTRODUCTION 


In the last decade, cosmology has passed into the cat¬ 
egory of precision sciences. Many cosmological parame¬ 
ters are currently determined with a high precision that 
occasionally reaches fractions of a percent (Ade et al. 
2014). One of such parameters is the baryon-to-photon 
ratio rj = nb/ti 7 , where rib and n 7 are the baryon and 
photon number densities in the Universe, respectively. 
In the standard cosmological model, the present value 
of 77 is assumed to have been formed upon completion of 
electron-positron annihilation several seconds after the 
Big Bang and has not changed up to now. 

The value of n 7 associated with the cosmic mi¬ 
crowave background (CMB) photons is defined by the 
well-known relation 
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where £(x) is the Riemann zeta function, k is the 
Boltzmann constant, h is the Planck constant, c is the 
speed of light, and T is the CMB temperature at the 
corresponding epoch. The CMB temperature is cur¬ 
rently determined with a high accuracy and is To = 
2.7255(6) K at the present epoch (Fixsen 2009); for other 
epochs, it is expressed by the relation T = T 0 (l + z), 
where z is the cosmological redshift at the corresponding 
epoch. Thus, given n 7 , a relation between the parameter 
77 and fib, the relative baryon density in the Universe, 
can be obtained (Steigman 2006): 

77 = 273.9 x 10" lo fl b /i 2 , 
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where h = 0.673(12) is the dimensionless Hubble param¬ 
eter at the present epoch (Ade et al. 2014). According 
to present views, the baryon density, which is the den¬ 
sity of ordinary matter (atoms, molecules, planets and 
stars, interstellar and intergalactic gases), does not ex¬ 
ceed 5% of the entire matter filling the Universe, while 
95% of the density in the Universe is composed of un¬ 
known forms of matter/energy that manifest themselves 
(for the time being) gravitationally (see, e.g., Gorbunov 
and Rubakov 2008). 

At present, observations allow fib to be indepen¬ 
dently estimated for four cosmological epochs: 

(i) the epoch of Big Bang nucleosynthesis (zbbn ~ 10 9 ; 
see, e.g., Steigman et al. 2007); 

(ii) the epoch of primordial recombination (zpr, — 1100; 
see, e.g., Ade et al. 2014); 

(iii) the epoch associated with the Lya forest (z ~ 24-3; 
i.e., ~10 Gyr ago; see, e.g., Rauch 1998; Hui et al. 2002); 

(iv) the present epoch (z = 0; see, e.g., Fukugita and 
Peebles 2004). 

For the processes at the epochs of Big Bang nucle¬ 
osynthesis and primordial recombination, r) is one of 
the key parameters determining their physics. For these 
epochs, the methods of estimating 77 , (i) comparing the 
observational data on the relative abundances of the 
primordial light elements (D, 4 He, 7 Li) with the pre¬ 
dictions of the Big Bang nucleosynthesis theory and (ii) 
analyzing the CMB anisotropy, give the most accurate 
estimates of 77 to date that coincide, within the observa¬ 
tional error limits: ?7bbn = (6.0 ±0.4) x 10 ~ 10 (Steigman 
2007) and ?7cmb = (6.05T0.07) x lO ' 10 (Ade et al. 2014). 
This argues for the correctness of the adopted model of 
the Universe and for the validity of the standard physics 
used in theoretical calculations. However, it should be 
noted that at present, as the accuracy of observations 
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increases, some discrepancy between the results of obser¬ 
vations and the abundances of the primordial elements 
predicted in the Big Bang nucleosynthesis theory has 
become evident. The “lithium problem” is well known 
(see, e.g., Cyburt et al. 2008); not all is ideal with he¬ 
lium and deuterium (for a detailed discussion of these 
problems, see Ivanchik et al. 2015). These inconsisten¬ 
cies can be related both to the systematic and statistical 
errors of experiments and to the manifestations of new 
physics (physics beyond the standard model). 

The determination of fib and the corresponding 77 
at epochs (iii) and (iv) has a considerably lower accu¬ 
racy. The value of 77 measured for the epoch associated 
with the Lya forest coincides, by an order of magnitude, 
with ? 7 bbn and 77 cmb> but, at the same time, is also 
strongly mo del-dependent (e.g., Hui et al. 2002). The 
measured 52b and 77 at the present epoch are at best 
half those predicted by Big Bang nucleosynthesis cal¬ 
culations and CMB anisotropy analysis. The so-called 
problem of missing baryons (see, e.g., Nicastro et al. 
2008) is associated with this. 

It is hoped that further observations and new exper¬ 
iments will allow fib for different cosmological epochs 
and the corresponding 77 to be determined with a higher 
accuracy. In turn, this can become a powerful tool for 
investigating the physics beyond the standard model, 
where the values of 77 for different cosmological epochs 
can be different. Constraints on the deviation of 77 allow 
various theoretical models admitting such a change to 
be selected. 

In this paper, we discuss the possibility of a change 
in 77 on cosmological time scales attributable to the de¬ 
cays of dark matter particles. For example, supersym¬ 
metric particles (see, e.g., Jungman et al. 1996; Bertone 
et al. 2004; and references therein) can act as such par¬ 
ticles; some of them can decay into the lightest stable 
supersymmetric particles and standard model particles 
(baryons, leptons, photons, etc.; see, e.g., Cirelli et al. 
2011 ): 


X^ X +- 


f 7 + 7 + ... 
\p + p + ..., 


(1) 


where X and \ are unstable and stable dark matter par¬ 
ticles, respectively. This can lead to a change in 77 . 

The currently available observational data suggest 
that the dark matter density in the Universe is approx¬ 
imately a factor of 5 larger than the baryon density: 
52 cdm — 552b, i.e., the relation between the number 
density of dark matter particles and the number densi¬ 
ties of baryons and photons in the Universe is 7 Tcdm — 
5(mb/77icDM)^b = 5(77Tb/7'ncDM)%f?- Assuming that 
the changes in the number densities of various types 
of particles in the decay reactions of dark matter parti¬ 
cles are related as Aucdm ~ A rib and A?7 cdm ~ An 7 , 
it is easy to see that the parameter 77 is most sen¬ 
sitive precisely to the change in baryon number den¬ 
sity. In the decays of dark matter particles with masses 
mcDM ~ lOGeV—ITeV, the change in 77 as a result 
of the change in baryon number density could reach 
A77/77 ~ 0.01 — Ifl The change in photon number density 

1 Here and below, out of all baryons, we restrict ourselves 
to protons. This assumption is valid for obtaining estimates, 
because the bulk of the baryon density in the Universe is 
contained in the hydrogen nuclei, while heavier baryons (for 


and the change in 77 attributable to it will be approxi¬ 
mately billion times smaller. Therefore, in our paper we 
focused our attention on the possibility of a change in 
?7 due to the decays of dark matter particles with the 
formation of a baryon component. 

Despite the negligible contribution to the change in 
77 from the photon component, a comparison of the pre¬ 
dicted gamma-ray background (dark matter particle de¬ 
cay products) with the observed isotropic gamma-ray 
background in the Universe can serve as an additional 
source of constraints on the decay models of dark mat¬ 
ter particles. The photons produced by such processes 
are high-energy ones. The observational data on the 
isotropic gamma-ray background constrain their possi¬ 
ble number in the Universe, which, in turn, narrows the 
range of admissible parameters of dark matter particles, 
determines the maximum possible number of baryons, 
the decay products of dark matter particles, and the cor¬ 
responding change in the baryon-to-photon ratio in such 
decays. Thus, the observational data on the gamma-ray 
background, along with the cosmological experiments 
described above, serve as a source of constraints on the 
decay models of dark matter particles and on the pos¬ 
sible change in 77 . Running ahead, we will say that at 
present the constraints from isotropic gamma-ray back¬ 
ground observations are more severe than those follow¬ 
ing from cosmological experiments. 

Depending on the lifetime of dark matter particles, a 
statistically significant change in 77 can occur at different 
cosmological epochs. We consider lifetimes r in the fol¬ 
lowing range: 2 bbn ^ t < to, where 2 bbn — 3 min is the 
age of the Universe at the end of the epoch of Big Bang 
nucleosynthesis, to — 13.8 Gyr is the present age of the 
Universe (Ade et al. 2014). The decays of dark matter 
particles with short lifetimes (r < 2 bbn) can change sig¬ 
nificantly the chemical composition of the Universe (see, 
e.g., Jedamzik 2004; Kawasaki et al. 2005). The available 
observational data on the abundances of the primordial 
light elements (D, 4 He, 7 Li) agree well with the predic¬ 
tions of Big Bang nucleosynthesis calculations, which, in 
turn, limits the possibility of such a change. For long life¬ 
times exceeding the present age of the Universe (r > to), 
the change in 77 at the above four cosmological epochs 
will be so small that this will unlikely allow it to be 
detected without a significant improvement in observa¬ 
tional capabilities. 

2. THE BARYON-TO-PHOTON RATIO IN 
MODELS WITH PARTICLE DECAY 

A large class of models with decaying dark matter par¬ 
ticles suggests the existence of the lightest stable parti¬ 
cle that we will designate as y. An unstable dark mat¬ 
ter particle, which we will designate as X, will decay 
with time into a y-particle and standard model par¬ 
ticles. There can be reactions of the type X —> ypp 
among such reactions, whose influence on 77 is investi¬ 
gated in this papei@. A quantitative parameter charac¬ 
terizing the fraction of the decay channels of X-particles 

example, D, He, etc.) are generated with a considerably lower 
probability. 

2 Since we consider cosmological time scales, all of the neu¬ 
trons and antineutrons that are also produced in such decays 
transform into protons and antiprotons. 
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whose products are hadrons (in our case, these will be 
protons and antiprotons) in the total number of decay 
channels is the hadronic branching ratio B} ll which is 
Bh = 1 in our case. 

The currently available observational data argue for 
the absence (or a negligible amount) of relic antimat¬ 
ter (baryon-asymmetric Universe). For this reason, the 
parameter 77 in the standard cosmological model is de¬ 
fined as the ratio of the baryon number density to the 
photon number density. Since in our model the decays 
of X-particles will lead to the production of protons and 
antiprotons, we will define the parameter 77 as the ra¬ 
tio of the sum of the baryon and antibaryon number 
densities to the photon number density in the Universe: 


3. INFLUENCE OF THE DECAY OF DARK 
MATTER PARTICLES ON THE CHANGE 
IN T) 

The evolution of the number densities of X-particles, 
%-particles, protons, and antiprotons in the Universe is 
described by the system of kinetic equations 

^^-+3 Hn x = — r n x , (5) 

at 

^-+3Hn x = Tn x , (6) 

nn — 

+ 3Hn PiP = —(av)pp n n p rip + B h Tn x , (7) 


V(z) = 


rih(z) + n B (z) 
n 7 (z) 


rib (z) + An p (z) + A n p (z) 


.BBN 


(z) 


( 2 ) 


?7bbn + A r](z), 


where ?r BBN and n BBN are the baryon and photon num¬ 
ber densities corresponding to 7 ?bbn = n BBN /n BBN ; 
A rip(z) and Arip(z) are the number densities of X- 
particle decay products: protons and antiprotons, re¬ 
spectively (in the model under consideration, An p (z) = 
A iip(z), i.e., the generated baryonic charge is A B = 0 ). 
It is this value of (2) that would be measured when 
determining the speed of sound of the baryon-photon 
plasma at the epoch of CMB anisotropy formation in the 
case of proton and antiproton generation in accordance 
with the formula (see, e.g., Gorbunov and Rubakov 
2010 ) 


2 = Sp = _ <?_ _ 

Sp 3 ( 1 + 3 p B g/ 4 p 7 ) 


( 3 ) 


where p B g = Pb + Pb is the sum of the baryon and an¬ 
tibaryon densities in the Universe. In the standard cos¬ 
mological model, this quantity coincides with the baryon 
density of the Universe p B - Thus, the baryon-to-photon 
ratio determined when analyzing the CMB anisotropy is 
also the ratio of the sum of the baryon and antibaryon 
number densities to the photon number density and has 
the following form in the presence of X-particle decay 
products: 


JJCMB = 


n h (z) +ng(z) 
n 7 (z) 


2 = ZPR 


??BBN + A? 7 (z P r), ( 4 ) 


Note that for very early decays the antiprotons be¬ 
ing produced have time to annihilate with protons, and 
77 again returns to its initial value 77 = ?7bbn- The de¬ 
cays of X-particles with long lifetimes will occur in an 
already fairly expanded Universe; consequently, the an¬ 
tiprotons being produced may not have time to annihi¬ 
late. Thus, the later 77 can differ from ?7bbn and ?7cmb- 
However, during the formation of a large-scale struc¬ 
ture, when halos in which the density of matter exceeds 
considerably the average one is formed, an excess of an¬ 
tiprotons would lead to enhanced gamma-ray radiation 
from them. 


where Eq. (7) consists of two equations describing the 
evolution of the proton and antiproton number densi¬ 
ties, ?i p and rip, respectively; nx and n x are the number 
densities of X- and %-particles, respectively; H = a/a is 
the Hubble parameter; a(t) is the scale factor; T = 1/r 
is the decay rate of X-particles; (crt')™ is the product 
of the relative velocity v and proton-antiproton annihi¬ 
lation cross section cr ann averaged over the momentum 
with a distribution function. In a wide energy range 
(lOMeV < T p < 10 GeV), this quantity may be con¬ 
sidered a constant, (tru) p p n = 10 - 15 cm 3 s -1 (see, e.g., 
Stecker 1967; Weniger et al. 2013). The parameters of 
the standard cosmological model presented in Table 1 
are used to solve Eqs. (5)-(7). 


Table 1. Cosmological parameters used in this paper 


Parameter 

Value 

Reference 1 

Hr 

5.46xl0 -5 

1 

ficDM 

0.265 

2 

fib 

0.05 

2 

Ha 

0.685 

2 

Ho 

67.3 kms -1 Mpc _1 

2 

to 

13.8 Gyr 

2 


1 1 - Fixsen (2009), 2 - Ade et al. (2014) 


Apart from the decays of dark matter particles, 
we investigated the processes of their annihilation. 
We showed that the influence of the annihilation of 
dark matter particles with an annihilation cross sec¬ 
tion (crn)™ = 10“ 26 cm 3 s _1 (see, e.g., Jungman et al. 
1996) in the case where the annihilation products are 
protons and antiprotons, XX PP, on the change in 
77 on all the time scales of interest could be neglected. 
This implies the absence of the terms responsible for 
the annihilation of X- and x-particles in Eq. (7). The 
change in 77 attributable to the annihilation of dark mat¬ 
ter particles with masses 10 GeV-1 TeV alone is negli¬ 
gible even at the epoch of Big Bang nucleosynthesis (at 
which the contribution from the annihilation is maxi¬ 
mal): |A?7/?7bbn| < 10~ 13 A 10 -11 (the upper limit cor¬ 
responds to a lower x-pai'ticle mass). 

To determine the initial conditions for Eqs. (5) and 
( 6 ), we introduce a parameter a defining the fraction 
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(by the number of particles) of unstable dark mat¬ 
ter particles in the entire dark matter at the epoch 
of Big Bang nucleosynthesis. For the range of lifetimes 
fuBN < t < to we consider, the entire dark matter at 
the present epoch will be composed of stable ^-particles 
some of which (ce) were produced by the decays of X- 
particles and some (1 — a) are the relic ones, i. e., the 
%-particle mass determines the initial conditions for the 
X-particles as well. The availability of reliable data on 
the parameter 77 at the epoch of Big Bang nucleosyn¬ 
thesis allows ??bbn to be used to determine the initial 
condition for Eq. (7). Thus, when solving the system of 
equations (5)-(7), we use the following initial conditions: 


0 „ -1 r »9 j .0 

z = ZBBN = 1U , t 


1 

2H(zbbn) ’ 


n 


0 

x 


rip = 77 bbn %( 2 : bbn ), rip = 0 , 


(!-«) 


^CDMPc 

2 5 
m x c z 


«x = a 


f^CDMPc 

2 ’ 
m x c z 


( 8 ) 


Let us write the system of equations (5)-(7) in a 
comoving volume that changes with time as ~ a 3 , i. e., 
~ (1 + z)~ 3 : 


dYx 

dt 

dYx 

dt 

dY p p 

dt 


—Tlx, 


ry x , 


-(av)^YpY p (l + z) 3 + B h TY x , 


(9) 

( 10 ) 

( 11 ) 


where Y) = rii/( 1 + z ) 3 is the number density of the itli 
type of particles in the comoving volume. 

In such a form, Eqs. (9) and (10) have obvious ana¬ 
lytical solutions that describe the evolution of the num¬ 
ber densities of X- and x-particles in the comoving vol¬ 
ume: 

Yx(t) = Yie-^ T , (12) 

Y x (t) = Y x ° + Y x °( l-e- t/T ), (13) 

where Y£ = n x /(l + 2 0 ) 3 and Y 3 = n°/( 1 + 2 0 ) 3 are 
the initial number densities of X- and x-particles in the 
comoving volume. Substituting solution (12), T = 1/r, 
and Bh = 1 into Eq. (11), we obtain the final system 
of equations describing the evolution of the proton and 
antiproton number densities in the model under consid¬ 
eration: 


dYp.p 

dt 


~{av) 


ann 

PP 


YpYpil + zf + ^-e-^. 


(14) 


The corresponding change in the baryon-to-photon 
ratio, 


A77(2) _ 77(2) - 7 ?BBN 


(15) 


??BBN ??BBN 


determined from the solution of the system of equations 
(14) for m x = lOGeV, a = 0.5, and various r is pre¬ 
sented in Fig. la. Note that the parameters a and m x 
enter into the system of equations (14) in the form of 
a ratio. Therefore, the result presented in Fig. 1 also 
corresponds to the case of larger masses of dark matter 
particles provided that a/m x is conserved. 

We see that the change in the baryon-to-photon ratio 
in the model under consideration for lifetimes r > 10 12 s 


can reach At7(2)/?7bbn ~ 0.01 — 1, which is a potentially 
observable value. We also see that the number densities 
of the protons and antiprotons in the comoving volume 
produced by late decays (r > 10 13 s) in an already fairly 
expanded Universe freeze in such a way that 77 can differ 
significantly from ?7bbn and 77 CMB by the present epoch. 
Note, however, that in the decays X —> ypp with the 
conservation of baryonic charge (i. e., An p (t) = Arip(f)), 
A?7/?7 bbn ~ 1 at the present epoch would imply al¬ 
most equal numbers of protons and antiprotons in the 
Universe, while our Universe is significantly asymmetric 
in baryonic charge. The existence of such a number of 
antiprotons in the Universe would also give rise to an 
excess of the gannna-ray background from the annihila¬ 
tion of protons with antiprotons (see the next section). 

Figure lb presents the dependence A?7(t)/?7bbn of 
the change in 77 at the epoch of primordial recombination 
(the epoch for which the parameter 77 has been measured 
most precisely to date) on the lifetime of X-particles r 
for various a. We see that the fraction of the change in 77 
at this epoch can reach A77/77BBN ~ 0.01 — 0.1, which is 
also a potentially observable value. Figure lc present the 
dependence A?7 (t)/77bbn referring to the present epoch 
(t 0 — 13.8 Gyr). We see that the decay of X-particles 
in the model under consideration leads to a significant 
change in the present baryon density for r > 10 13 s. 
However, the accuracy of its determination at an epoch 
2 ~ 2 — 3 and at the present epoch is still considerably 
lower than that for the epochs of Big Bang nucleosyn¬ 
thesis and primordial recombination. 

The results obtained should not come into conflict 
with other observational data: 

(1) The decays with a predominance of hadronic chan¬ 
nels at early epochs r <C <pr can change significantly 
the chemical composition of the Universe (see, e.g., 
Jedamzik 2004; Kawasaki et al. 2005). The available ob¬ 
servational data on the abundances of the primordial 
light elements (D, 4 He, 7 Li) agree well with the predic¬ 
tions of Big Bang nucleosynthesis calculations, which, 
in turn, limits the possibility of such a change. 

(2) The decays with r ~ tpR can distort the CMB spec¬ 
trum and affect the angular CMB anisotropy (see, e.g., 
Chen and Kamionkowski 2004; Chluba and Sunyaev, 
2012). Comparison with observational data also allows 
the possible models to be constrained severely. 

(3) The hadronic decays with r > tpR can give rise to 
an excess gamma-ray background from the annihilation 
of produced antiprotons with background protons and 
directly from the decays of X-particles (see the next sec¬ 
tion). 

In our case, we used data on the isotropic gamma- 
ray background to obtain constraints on the decays of 
particles with tpR < r < to, because a maximal effect 
of change in the baryon-to-photon ratio is expected for 
such lifetimes of X-particles (see Fig. 1). As we will see, 
at present these constraints are more significant than 
those that can be given by present-day cosmological ex¬ 
periments. 
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10" 9 10' 8 10" 7 10' 6 10' 5 10" 4 10' 3 10' 2 10' 1 10° 

1 /(1 + z) 



10 11 10 12 10 13 10 14 10 15 10 16 10 17 


X, s 

Figure 1. Fraction of the change in the baryon-to-photon ratio attributable to the decays of 

X-particles with lifetimes 10 5 s < r < 10 17 s (m x = lOGeV, a = 0.5); the vertical lines mark the epochs of Big 
Bang nucleosynthesis (zbbn ~ 10 9 ) and primordial recombination (zpr ~ 1100). The dependence At](t)/ijbbn of 
the change in the baryon-to-photon ratio at the epoch of primordial recombination t = tpR (b) and at the present 
epoch t = to (c) on the lifetime of X-particles for various values of the parameter a. 
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Figure 2. Isotropic gamma-ray background d<^ 1 /dVLdE 1 attributable directly to the decays of dark matter particles 
(dotted curve) and the annihilation of protons with antiprotons (dash-dotted curve) in the decay model of X- 
particles with a lifetime r = 10 14 s under consideration (m x = 10 GeV, mx — m x = lOGeV, a = 5 • 10 -6 ). The 
squares, circles, and triangles mark the experimental data taken from Ajello et al. (2008), Bloemen et al. (1999), 
and Sreekumar et al. (1998), respectively; the solid curve represents a fit to the experimental data from Gruber 
et al. (1999). 


4. CONSTRAINT ON THE POSSIBLE 

CHANGE IN T] ASSOCIATED WITH THE 
OBSERVATION OF AN ISOTROPIC 
GAMMA-RAY BACKGROUND 

As was shown by Cirelli et al. (2011), apart from pro¬ 
tons and antiprotons, photons and leptons will also be 
present among the end decay products of dark matter 
particles, with their fraction exceeding considerably the 
fraction of baryons even in the case of Bh = 1 (i. e., 
when the decays completely run via hadronic channels). 
The reason is that apart from protons and antipro¬ 
tons, mesons are produced in the hadronization process, 
which contribute to the photon and lepton components. 
In addition, the appearance of an antiproton fraction 
in the Universe will be accompanied by the formation 
of an additional gamma-ray background from the anni¬ 
hilation of proton-antiproton pairs. The main gamma- 
ray background as a result of such a process will arise 
from the decay of the 7r° meson produced by the proton- 
antiproton annihilation (Stecker 1967; Steigman 1976). 
Both these processes, which can be represented schemat¬ 
ically as 


{ 7 + 7+ ••• 

J^° -+ 7 + 7 ( 16 ) 

P P ->■ ^ ± + ^m( i +)> 

/J ± e 4 + t'e(z'e) + *+(*+)’ 

will contribute to the isotropic gamma-ray background 
in the Universe. 

We calculate the corresponding gamma-ray back¬ 
ground by taking into account its extension to cosmo¬ 
logical distances. Note that photons of different ener¬ 
gies at different cosmological epochs interact differently 
with the medium in which they propagate (see, e. g., 


Zdziarski and Svensson 1989; Chen and Kamionkowski 
2004). More specifically, there is a transparency win¬ 
dow: the photons with energies U 7 < 10 GeV emitted at 
epochs 0 < z < 1000 propagate almost without absorp¬ 
tion and reach us in the form of an isotropic gamma-ray 
background. The formation of such a gamma-ray back¬ 
ground is expected from the decays of X-particles with 
lifetimes tpR < r < to. 

The general formula describing the intensity 
of the isotropic gamma-ray background J 7 (U 7 ) 
(keV-cm _2 s~ 1 sr _1 keV _1 ) from various processes is 
(see, e.g., Peacock 2010) 


I-y(Ej) = Ej 


ci<I » 7 

dttdEj 


(17) 


c 

47 T 


1000 


dz- 


y([l + z \E~i) z ) c - t (e^,z) 


H(z)( 1 


z) 4 


where <f> 7 is the gamma-ray photon flux per unit time 
through a unit area, t{E 1 ,z ) is the optical depth de¬ 
scribing the absorption of a photon emitted at epoch 2 
with energy E 7 (l + z), e 7 is the volume emissivity that 
in our case is the sum of two terms, 

e 7 (U 7 , z) = e*(U 7 , z) + ePP(U 7 , z), (18) 

describing the two contributions to the gamma-ray 
background mentioned above. The first term e 7 is re¬ 
lated to the photons that are the X-particle decay prod¬ 
ucts; the second term is related to the photons that 
are the proton-antiproton annihilation products. These 
terms are described by the expressions 

^j{E 1 ,z) = E 1 Tn x (z)^ L (19) 

= E 7 TY x (z)(1 + z) 3 ^-, 
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Figure 3. Isotropic gamma-ray background d<& 1 / d£ldE 1 in the decay model of X-particles with lifetimes 10 14 s < 
t < 10 17 s under consideration for a = 1 (dash-dotted curve), a = 10 -3 (dotted curve), and a = a max (dashed 
curve) ( m x = 10 GeV, mx — rn x = 10 GeV). The solid curve represents a fit to the experimental data from Gruber 
et al. (1999). 


Table 2. Maximum admissible fraction of X-particles a max with various lifetimes r for x-partide masses of 10, 
100, and 1000 GeV and the corresponding maximum admissible change in the baryon-to-photon ratio A^/t/bbn 
at epoch z*. 


r, s 

C^max 

Ar/(z*) 

?7BBN 

z* 

m x = 10 GeV 

m x = 100 GeV 

m x = 1000 GeV 

10 14 

5 x 10 -6 

5 x 10" 5 

5 x 10" 4 

2.3 x 10 -6 

120 

10 15 

5 x 10" 7 

5 x 10 -6 

5 x 10 -5 

4.2 x 10 -7 

18 

10 16 

10' 7 

10“ 6 

10 -5 

10- 7 

2.2 

10 17 

10~ 8 

10- 7 

10“ 6 

10 -8 

0 


tf P ( E ~n z) = E i(°v) a ™n v (z)n p (z)-^ ( 20 ) 

= E y (av) p f*Y p (z)Y p (z)(l + z) 6 ^, 

where dN 1 /dE 1 is the spectrum of the photons 
(phot-keV _1 ) emitted in one event of X-particle de¬ 
cay (in Eq. (19)) and proton-antiproton annihilation (in 
Eq. (20)). 

In our calculations, we use the spectrum dN x /dE x of 
the photons that are the X-particle decay products cal¬ 
culated in the PYTHIA package. The numerical code 
for computing the spectra of the dark matter particle 
decay and annihilation products was taken from the 
sitqj; the details of using it can be found in Cirelli et 

3 http://www.marcocirelli.net/PPPC4DMID.html 


al. (2011). We use the data for tox — m x ~ 10 GeV 
from the entire range of energy release accessible in the 
numerical code in such reactions, 10 GeV—200 TeV, to 
determine an upper bound on the possible change in 
77 . The optical depths in (17) were also taken from this 
site. The spectrum dN 1 /dE 1 of the photons that are 
the proton-antiproton annihilation products was taken 
from Backenstoss et al. (1983). 

For comparison, Fig. 2 presents the gamma-ray 
background d<b 7 / dfldE^ (phot-cm _ 2 s _ 1 sr _ 1 keV _1 ) at¬ 
tributable to the contribution from each of the two terms 
in (18). The observational data on the isotropic gamma- 
ray background (lOkeV—1 GeV) taken from Sreekumar 
et al. (1998), Bloemen et al. (1999), Gruber et al. (1999), 
and Ajcllo et al. (2008) are also presented in the figure. 
We see that the gamma-ray background directly from 
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Figure 4. Fraction of the change in the baryon-to-photon ratio A ' l<z> = — 222 tL attributable to the decays of 
X-particles with lifetimes 10 14 s < t < 10 17 s (m x = lOGeV, TOx — m x = lOGeV). The solid and dashed curves 
correspond to the fractions of X-particles a = 1 and a = a max , respectively. The vertical line marks the epoch of 
primordial recombination (zpr, — 1100 ). 


the decays of X-particles allow stringent constraints to 
be placed on the decay processes. 

Figure 3 shows the total gamma-ray background 
d<f > 7 /d£ldE 1 with the inclusion of both terms in (18) 
for X-particle lifetimes fpR < r < to and various val¬ 
ues of the parameter a. The gamma-ray background 
admissible by the currently available observational data 
corresponds to a max , which characterizes the maximum 
admissible fraction of unstable X-particles with the cor¬ 
responding lifetime. The values of a max for lifetimes 
tpR r < to are presented in Table 2. 

Figure 4 presents the fraction of the change in the 
baryon-to-photon ratio corresponding to a max for var¬ 
ious X-particle lifetimes (for comparison, Fig. 4 also 
presents this change for a = 1). We see that this change 
may reach Ati(z)/iibbn < 10 -5 . The present-day obser¬ 
vational accuracy is A 77/77 ~ 1CD 2 —10 _1 . Note that the 
corresponding number of antiprotons in the Universe at 
the present epoch related to A?; via the relation 


np 

n v 


1 Ary 

2 r?BBN 


2=0 


is consistent with the observational data on antiprotons 
in cosmic rays (see, e.g., Adriani et al. 2010). 

Since the parameters a and m x enter into the sys¬ 
tem of equations (14) in the form of a ratio, the result 
obtained can be easily generalized to the case of larger 
masses of dark matter particles. For x-particles with 
masses m x = 10, 100, and 1000 GeV, the derived pa¬ 
rameter a max and the corresponding maximum change 
Ajy/ryBBN in the baryon-to-photon ratio are listed in 


Table 2. The table also gives the cosmological redshift 
z* corresponding to the maximum change in 77. 


5. CONCLUSIONS 

We investigated the influence of the baryonic decay 
channels of dark matter particles X —> ypp on the 
change in the baryon-to-photon ratio at different cos¬ 
mological epochs. 

We showed that the present dark matter density 
ficDM — 0.26 is sufficient for the decay reactions of dark 
matter particles with masses lOGeV—ITeV to change 
the baryon-to-photon ration up to At 7 (z)/? 7 bbn ~ 0.01 — 
1 (Fig 1). However, such a change in 77 would lead to an 
excess of the gannna-ray background from the annihi¬ 
lation of proton-antiproton pairs, the decay products of 
dark matter particles, and from the gamma-ray photons 
produced directly in the decays of dark matter particles. 

We used the observational data on the isotropic 
gamma-ray background to constrain the decay models 
of dark matter particles leading to a maximum effect 
of change in 77 : we determined the maximum admissible 
fraction of unstable dark matter particles with lifetimes 
fpR < r < to and the change in 77 related to them. The 
maximum possible change in the baryon-to-photon ra¬ 
tio attributable to such decays is Ar](z)/r)BBN 10 -5 
(Fig. 4). 

Despite the fact that at present the data on the 
gamma-ray background constrain most severely the de¬ 
cay models of dark matter particles with the emission 
of baryons, the situation can change in future, with in¬ 
creasing accuracy of existing cosmological experiments 
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and the appearance of new ones. The detection of a 
change in the baryon-to-photon ratio in such experi¬ 
ments at a level of < 10 -5 will serve as evidence for 
the existence of decaying dark matter particles, while its 
detailed study will be a powerful tool for studying their 
properties. In contrast, the constancy of the baryon-to- 
photon ratio will serve as a new source of constraints 
on the range of admissible parameters of dark matter 
particles. 
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